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Resumen — 

En este trabajo proponemos una VNS eficiente pa¬ 
ra el problema del bosque generador de /c-etiquetado. 
Este problema es una extensión del problema del 
árbol generador de mínimo etiquetado ( Mínimum La- 
belling Spanning Tree Problem) con importantes aplica¬ 
ciones en redes de telecomunicaciones y transporte 
multimodal. Se trata de, dado un grafo no dirigido 
cuyos enlaces están etiquetados y un número ente¬ 
ro positivo k encontrar el bosque o subgrafo con el 
menor número de componentes conexas usando a lo 
sumo k etiquetas distintas. Para abordar lo propone 
una VNS básica (BVNS, Basic Variable Neighbourhood 
Search) donde la amplitud máxima del proceso de agi¬ 
tación es el parámetro clave. Se estudian diferentes es¬ 
trategias para establecer el valor de dicho parámetro 
obteniéndose como mejor estrategia hacerlo depender 
del estado del proceso de búsqueda. La BVNS con la 
estrategia seleccionada es experimentalmente compa¬ 
rada con las metaheurísticas que han aparecido en la 
literatura aplicada a este tipo de problema. 

Palabras clave —VNS, Metaheurísticas, Grafos, Bos¬ 
que Generador de /c-etiquetado. 


I. Introducción 


blema del bosque generador fc-etiquetado. Se denota 
por fcLSFP (fc-¿abelled Spanning Forest Problem) y 
se formaliza en los siguientes términos: 

Sea fc un entero positivo y G = (V. A, L) un gra¬ 
fo no dirigido etiquetado donde V es el conjunto 
de sus vértices y A el de sus aristas que están eti¬ 
quetadas con etiquetas del conjunto finito de eti¬ 
quetas L. Se trata de encontrar el conjunto de eti¬ 
quetas L* C L de no más de k etiquetas distin¬ 
tas tal que el subgrafo G* de G tomando todas las 
aristas cuya etiquetas están en L* tenga el menor 
número posible de componentes conexas. Formal¬ 
mente, denotamos por l(a) € L la etiqueta de la 
arista a € A. Para cualquier conjunto de etiquetas 
C C L denotemos por G(C) al subgrafo del gra¬ 
fo etiquetado G formado por todas las aristas con 
etiquetas en C; es decir G{C) = (V,A(C)) don¬ 
de A{C) = {a e A : l(a) G C}. Denotemos por 
Comp(C) el número de componentes conexas del 
subgrafo G(C). El problema es: 


Las redes de transporte multimodal pueden repre¬ 
sentarse por grafos donde cada enlace tiene una eti¬ 
queta que denota una empresa que lo gestiona o ti¬ 
po de transporte diferente. En muchos contextos de 
transporte multimodal se pretende garantizar la ma¬ 
yor conectividad posible entre los nodos de la red uti¬ 
lizando para ello el número mínimo de empresas pro¬ 
veedoras o tipos de transporte El, El- El propósi¬ 
to es reducir el coste de la implantación del siste¬ 
ma de transporte o la complejidad global de la red. 
En algunas situaciones reales estaremos interesados 
en la optimización de la conectividad con un límite 
en el número de proveedores o modos de transporte 
utilizados 1141 . Este problema denominado problema 
del bosque generador etiquetado (Labelled Spanning 
Forest Problem, LSFP), propuesto originalmente en 
[2], pertenece a una clase de problemas derivados del 
problema del árbol generador de etiquetado míni¬ 
mo ( Mínimum Labelling Spanning Tree Problem , 
MLSTP) [5], que ha recibido bastante atención re¬ 
cientemente. El problema del bosque generador eti¬ 
quetado, con un entero positivo k como cota superior 
del número de etiquetas a utilizar, se denomina pro- 
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mín{Comp{C) : C C L, \C\ < k}. 

Como solución al problema se aporta, junto al con¬ 
junto de etiquetas óptimo L*, el grafo asociado G* = 
G(V,A(L*)) o un bosque generador H*\ es decir un 
subgrafo con el mismo número de componentes co¬ 
nexas que G* pero sin ciclos. Si la solución óptima 
G* contiene ciclos se pueden romper arbitrariamente 
cada uno de ellos quitando una arista cualquiera en 
tiempo polinomial para obtener el bosque generador 
H*. 

En [5] se prueba que el problema del bosque ge¬ 
nerador /c-etiquetado, /cLSFP, es un problema NP- 
duro. El fcLSFP generaliza el problema del árbol ge¬ 
nerador de mínimo etiquetado (MLSTP, Mínimum 
Labelling Spanning Tree Problem ) cuyo objetivo es 
determinar el árbol generador del grafo no dirigi¬ 
do etiquetado con el número mínimo de etiquetas. 
Una solución al problema MLST sería también una 
solución al problema fcLSF si el árbol solución ob¬ 
tenido no sobrepasa el límite de fc en el número de 
etiquetas utilizadas. Por lo tanto los métodos de so¬ 
lución del problema fcLSF se deben buscar en pri¬ 
mer lugar adaptando los propuestos para el MLSTP. 
En ¡5] se propuso el algoritmo MVCA ( Máximum 
Vertex Covering Algorithm ) como primera heurística 
para el problema MLST que fue mejorada por El; 
otras propuestas metaheurísticas para dicho proble¬ 
ma aparecen en mi, a, 0,0,0, m , i, m, 

ra- 


La estructura del resto del trabajo es como si¬ 
gue. En la Sección [Ti] se presentan los detalles de 
los métodos de solución para el fcLSFP aparecidos 
en la literatura propuestos por diversos autores. En 
la Sección [ITT] se describe con algo más de detalle la 
VNS básica propuesta y el ajuste de la amplitud de 
agitación para el problema. La Sección IIVI muestra 
algunos resultados de la comparación experimental 
y finalizamos con la sección de conclusiones fVl 

II. Algoritmos para el fcLSFP 

Se describe, en primer lugar, un método exacto 
derivado de [7 para el kLSF. A continuación se 
describen las principales metalieurísticas aparecidas 
en la literatura: un Método Piloto 0 , m y un algo¬ 
ritmo genético Se presenta la adaptación al 

fcLSFP del procedimiento de búsqueda voraz aleato- 
rizada adaptativa GRASP propuesta en [7] para el 
MLSTP. 

Todos estos algoritmos se basan, de una mane¬ 
ra u otra, en el algoritmo MVCA 0 , m adapta¬ 
do al fcLSF por 2] como subprocedimiento cons¬ 
tructivo. A partir de un conjunto inicialmente vacío 
de etiquetas C, el MVCA añade iterativamente a 
C la etiqueta l £ (L — C) tal que el subgrafo 
G(C) = (V,A(C)) tenga el número mínimo de com¬ 
ponentes conexas \Comp(C)\, y se detiene cuando el 
número de etiquetas utilizadas alcance la cota supe¬ 
rior (|Cj = fc), o cuando quede sólo una componente 
conexa {\Comp{C)\ = 1). 

A. Método Exacto 

El enfoque exacto para el fcLSF se basa en un pro¬ 
ceso de exploración por vuelta atrás (backtracking). 
Dado el grafo no dirigido etiquetado G = (V,A,L) 
con n vértices, m aristas y i etiquetas, y un número 
entero positivo fc como límite superior para el núme¬ 
ro de etiquetas, el método exacto realiza una rami¬ 
ficación y poda en el espacio de soluciones parciales 
basada en un procedimiento recursivo. Se parte de 
un conjunto de etiquetas vacío y se construye una 
solución mediante la adición iterativa de etiquetas 
hasta que el número de etiquetas llegue a fc. Cuando 
se alcanza la cota se aplica la fase de vuelta atrás 
para explorar nuevas alternativas. El algoritmo se 
detiene si se encuentra una solución con sólo una 
componente conexa sin superar el número fc de eti¬ 
quetas utilizadas. Este proceso da lugar a una enu¬ 
meración de todas las posibles combinaciones de fc 
etiquetas. Por lo tanto, la combinación con el menor 
número de componentes conexas es la solución exac¬ 
ta proporcionada como resultado del método. Aun¬ 
que el tiempo de cómputo de este procedimiento es 
exponencial en el número de etiquetas, si el tamaño 
del problema es moderado o el valor de la solución 
óptima es pequeño, el tiempo de ejecución de este 
método exacto es razonable y es posible obtener la 
solución exacta. 


B. Método Piloto 

La idea central del Método Piloto (PM) 0, [2] 
consiste en agotar todas las opciones posibles con 
respecto a una llamada solución maestra, por me¬ 
dio de una heurística básica. Esta heurística básica 
realiza tentativamente iteraciones con respecto a la 
solución maestra hasta que se evalúan todas las po¬ 
sibles elecciones locales. La mejor solución hasta en¬ 
tonces se convierte en la nueva solución maestra, y 
el procedimiento continúa hasta que se cumplan las 
condiciones de parada predefinidas, como alcanzar 
cierto tiempo máximo permitido de CPU o cierto 
número máximo de iteraciones. 

El Método Piloto para el problema fcLSF parte 
de la solución nula (un conjunto vacío de etiquetas) 
como solución maestra M. Entonces, para cada eti¬ 
queta i £ M, se trata tentativamente de extender 
una copia de M a una solución completa incluyendo 
la etiqueta i , construida mediante el uso del algo¬ 
ritmo MVCA como heurística básica. Se inserta en 
la solución parcial la etiqueta que más disminuye el 
número componentes conexas del subgrafo parcial 
con una estrategia voraz ( greedy ) ¡2j. El número de 
componentes conexas de la solución obtenida a par¬ 
tir de M <-MU{j} se usa como función objetivo de 
cada candidato i ^ M. Cuando todas las etiquetas 
candidatas posibles con respecto a la solución maes¬ 
tra han sido evaluadas, la etiqueta candidata i* con 
menor valor de la función objetivo se añade a la solu¬ 
ción maestra (M «— M U {**}). Sobre la base de esta 
nueva solución maestra Ai, se inician nuevas itera¬ 
ciones Vi ^ Ai, proporcionando un nuevo elemento 
de la solución de i*, y así sucesivamente. Este meca¬ 
nismo se itera hasta que no se necesite añadir más 
etiquetas a la solución maestra (es decir, hasta que 
el número de etiquetas utilizadas |AI| alcanza fc, o se 
obtenga una única componente conexa). Alternati¬ 
vamente, se pueden imponer algunas condiciones de 
finalización que permitan que el algoritmo pueda al¬ 
canzar estas condiciones. La última solución maestra 
corresponde a la mejor solución hasta entonces y se 
aporta como respuesta del método. 

C. Algoritmo Genético 

Los Algoritmos Genéticos (AGs) se inspiran en 
los principios de la selección natural, la adaptación 
al entorno, la evolución y operaciones de cruce y 
mutación entre los individuos dentro de una pobla¬ 
ción [roj. Para la aplicación de esta metaheurísti- 
ca al fcLSFP, dado el grafo etiquetado no dirigido 
G = ( V,A,L ), un individuo (o cromosoma) de la 
población es un subgrafo de G, donde cada posible 
etiqueta se corresponde un gen, y el grado de adap¬ 
tación ( fitness ) se define como el número de compo¬ 
nentes conexas ¡mi- La población inicial se gene¬ 
ra añadiendo etiquetas al azar a un conjunto vacío, 
hasta que aparecen individuos con fc etiquetas. Los 
operadores de cruce y mutación se aplican a conti¬ 
nuación para dar lugar a una nueva generación desde 


la anterior m, m ■ Las probabilidades de cruce y 
mutación se establecen al 100 % m, 0- El número 
de generaciones se establece como la mitad del ta¬ 
maño de la población. Por lo tanto, en este algoritmo 
genético cu, m el único parámetro a ajustar es el 
número de generaciones o el tamaño de la población. 

El cruzamiento da lugar a un único descendiente 
para cada dos soluciones factibles padres. Dados los 
dos padres P\ C L y P 2 C L, se comienza formando 
su unión P = P\ U P 2 . A continuación se agregan 
etiquetas de P a una solución inicialmente vacía uti¬ 
lizando la estrategia del MVCA hasta que se obtenga 
una solución hija con fc etiquetas. Es decir, se aplica 
el MVCA al subgrafo G(P ) de etiquetas en P. Por 
otro lado, la operación de mutación consiste en una 
perturbación: se añade una nueva etiqueta a al azar, 
y se elimina la etiqueta que dé lugar al individuo con 
mejor ajuste. 

Tras de una serie de generaciones el algoritmo con¬ 
verge y el mejor individuo de la población se pro¬ 
porciona como resultado. Téngase en cuenta que el 
algoritmo genético termina antes de la última gene¬ 
ración si se encuentra una solución factible formada 
por una sola componente conexa. 

D. Búsqueda Adaptativa Voraz Aletaorizada 

El procedimiento de búsqueda adaptativa voraz 
aletaorizada conocida como GRASP (Greedy Ran- 
domized Adaptiva Search Procedure) es básicamente 
una metaheurística constructiva que consta de una 
fase de iterativa constructiva seguida de una fase de 
mejora basada en una búsqueda local. En la fase de 
constructiva se construye una solución mediante la 
aplicación de un procedimiento voraz aleatorizado. 
Se crea iterativamente una lista restringida de candi¬ 
datos formada por los candidatos más prometedores 
para ser incorporados a la solución parcial inicial¬ 
mente vacía. La lista restringida de candidatos se 
denota por LRC y, en cada iteración, se selecciona 
de ella al azar uno de sus elementos que se incorpo¬ 
ra a la solución parcial [15] • Para el fcLSFP utiliza¬ 
mos de una lista restringida de candidatos basada 
en la evaluación de las etiquetas que se incorporan a 
LRC a [7]. Es una extensión del criterio clásico apli¬ 
cado en GRASP, que consiste en colocar en la lista 
sólo las etiquetas candidatas que tienen una valora¬ 
ción (el número de componentes conexas en el caso 
del fcLSFP) no mayor que un umbral definido por el 
usuario m- 

En nuestra implementación, se aplica una alea- 
torización completa para elegir la etiqueta inicial a 
incorporar la solución vacía. Esto corresponde a es¬ 
tablecer en la primera iteración el umbral en + 00 , lo 
que significa que la lista de candidatos está formada 
por todas las etiquetas del modelo. Para las restan¬ 
tes iteraciones, para cada nueva etiqueta a añadir, la 
lista LRC se forma considerando sólo las etiquetas 
que dan el número mínimo de componentes conexas 
en la etapa específica, con el fin de intensificar al 


máximo el proceso de búsqueda. Esto significa que 
se fija el umbral al número mínimo de componen¬ 
tes conexas obtenidas por las posibles etiquetas a 
añadir [7]. Es decir, las etiquetas que dan lugar al 
menor número de componentes conexas en cada paso 
constituyen la lista restringida de candidatos LRC. 
Luego, en cada iteración se selecciona aleatoriamen¬ 
te una etiqueta de la lista restringida de candidatos 
LRC que se añade a la solución actual. Se actuali¬ 
za la lista restringida de candidatos para la nueva 
solución parcial y se itera el proceso. La fase cons¬ 
tructiva finaliza cuando se obtenga una solución con 
fc etiquetas, o que dé lugar a un subgrafo con solo 
una componente conexa. 

La consiguiente mejora por búsqueda local quita 
con criterio greedy una etiqueta y añade otra etique¬ 
ta utilizando el algoritmo MVCA, de forma que el 
número total de componentes conexas se reduzca al 
mínimo. Esto produce una mayor intensificación del 
algoritmo. 

Este proceso de dos fases es iterativo, continuando 
hasta que se satisfaga la condición de parada especi¬ 
ficada por el usuario. Los criterios de parada usuales 
hacen referencia a que se alcance un tiempo de CPU 
máximo permitido, a un número máximo de itera¬ 
ciones, o a un número máximo de iteraciones entre 
dos mejoras sucesivas. El resultado final aportado 
por el procedimiento GRASP es la mejor solución 
encontrada hasta entonces en todas las iteraciones. 

III. La Búsqueda por Entornos Variables 

La idea fundamental en una Búsqueda por 
Entorno Variable (VNS, Variable Neighbourhood 
Search ) es el cambio sistemático en la estructura de 
entornos para escapar de los mínimos locales m- 
La idea de la Búsqueda por Entorno Variable Bási¬ 
ca (BVNS, Basic Variable Neighbourhood Search ) es 
explorar entornos cada vez más alejados para iniciar 
la búsqueda local y así tratar de mejorar el ópti¬ 
mo local alcanzado ¡12j . Nuestra implementación 
de BVNS para el fcLSFP está inspirada en el éxi¬ 
to de la VNS propuesta para el problema de árbol 
generador de etiquetado mínimo MLST ( Mínimum 
Labeling Spanning Tree ) [7]. Dado el grafo no di¬ 

rigido G = (V, A) con \V\ = n vértices y \A\ = m 
aristas que están etiquetadas con \L\ = í etiquetas 
y un entero positivo fc como límite para el número 
de etiquetas a utilizar, cada solución C C L se co¬ 
difica por una cadena binaria ( String ) de la forma 
C = [h,l 2 ,..., I eii], donde U = 1 si la etiqueta i se in¬ 
cluye en la solución C, y = 0 en caso contrario. Se 
selecciona el conjunto de q m ax estructuras de entor¬ 
nos N q , con g=l,2,..., q m axi cada una de las cuales 
es una aplicación del espacio de soluciones en las par¬ 
tes o subconjuntos de soluciones. El parámetro q m ax 
es el número máximo de estructuras de entornos a 
utilizar. 

La elección más común y más simple de estruturas 
de entornos para esta metaheurística es una fami- 


lia de entornos de tamaño creciente; es decir, tales 
que | -/Vi ((7)1 < ... < \N q (C)\ < ... < \N qmax (C)l 
VC C L. Sea p(C 1: C 2 ) = ¡CiACal = sum¿[ A¿ la 
distancia de Hamming entre dos soluciones C\ y C 2 , 
donde A¿ = 1 si la etiqueta i se incluye en una de 
las dos soluciones, pero no en la otra, y A¡ = 0 en 
caso contrario, Vi = 1,2, El conjunto de solu¬ 
ciones del g-ésimo entorno de una solución de C se 
establece por N q (C) = {C C L : p(C',C) = q}, 

Vq = 1, ...,5mai- 

La VNS básica que hemos implementado comienza 
por una solución inicial C*. Este conjunto inicial de 
etiquetas se obtiene mediante la inclusión iterativa 
de etiquetas seleccionadas de forma aleatoria en un 
conjunto inicialmente vacío hasta tener k etiquetas. 
A continuación, el algoritmo aplica la fase de agita¬ 
ción desde C* , que consiste en la selección aleatoria 
de una solución C' del entorno N q (C*) de la solu¬ 
ción actual C*, dejando variar al parámetro q desde 
1 hasta qmax a lo largo de la ejecución. Para construir 
una solución C' £ N q (C*), el proceso de agitación 
comienza con C' = C*. En primer lugar, se quitan q 
etiquetas al azar de C", y, si después de que las eti¬ 
quetas hayan sido retiradas, se requiere una mayor 
expansión de la estructura de entorno (es decir, en 
el caso de que q > |C"|) se incluyen nuevas etique¬ 
tas adicionales en C' elegidas al azar del conjunto de 
etiquetas sin usar ( L — C'), [7]. El pseudocódigo del 
proceso de agitación se describe en el Algoritmo [Tj 


se alcance k etiquetas. Con esta búsqueda local, 
se obtiene un nuevo conjunto C" de etiquetas. Si 
con C" no se ha obtenido una mejora (es decir, si 
\Comp(C")\ > |Co77Tp(C*)|, la amplitud de la agi¬ 
tación se incrementa (q <— q + 1), dando lugar a 
una mayor diversificación del proceso de búsqueda. 
El pseudocódigo de la búsqueda local se describe en 
el Algoritmo [2j 


Algorithm 2: Búsqueda local de la BVNS para el 
fcLSFP._ 

Procedimiento Búsqueda-local(C"): 

begin 

Sea C" «- C"; 
i t— 1; 

while i < \C'\ do 

[Eliminar la etiqueta í-ésima de C"; sea C": 

C" <-C- {c¿}; 

while ( Comp(C ") > 1) & (|Cj < k ) do 
Seleccionar al azar una etiqueta no usada 
c £ (L — C") que minimice Comp(C" U {c}); 
Añadir la etiqueta c al conjunto de etiquetas 
usadas: C" t— C" U {c}; 

Actualizar H" = (V,A{C")) y Comp(C "); 
end 

i £- i + 1 ; 

end 

end 


Algorithm 1: Agitación de la BVNS para el fcLSFP. 
Procedimiento Agitación(iV g (C')): 

begin 

Sea C' £~ C\ 

1 ^— 1 ; 

while i < q do 
if i < \C\ then 

Seleccionar al azar una etiqueta c' £ C'\ 
Borrar la etiqueta c' del conjunto de etiquetas 
usadas: C £- C — {d} ; 
else 

Seleccionar al azar una etiqueta no usada c', i. 
d £ {L - C ); 

Añadir la etiqueta d al conjunto de etiquetas 
usadas: C' <— C’ U {c'}; 

end 

z ^— z -f- 1 ; 

end 

Actualizar H’ = (V, A{C')) y Comp(C')\ 

end 


La VNS básica (BVNS T2] comienza con una so¬ 
lución inicial C*. Se aplica el procedimiento de agi¬ 
tación descrito en el Algoritmo Q] desde C* para ob¬ 
tener C’ £ N q (C*) y desde C' como solución inicial 
se ejecuta la búsqueda local descrita en Algoritmo 
[2] para obtener C". Si C" mejora a C* (es decir, si 
| Comp{C")\ < |Comp(C'*)|) la sustituye como solu¬ 
ción actual y se reinicia la amplitud q de la agitación 
a 1; (q t— 1). En otro caso, se incrementa en una 
unidad la amplitud de la agitación; q •<— q + 1. De 
esta forma la amplitud de la agitación q varía desde 
1 hasta un valor máximo q m ax que es el parámetro 
eclave de la BVNS. El proceso de cambiar la estructu¬ 
ra de entornos cuando la búsqueda local alcanza un 
mínimo local que no mejora la mejor solución hasta 
entonces representa la idea central de la VNS bási¬ 
ca al proporcionar una diversificación creciente. Este 
procedimiento se repite iterativamente hasta que se 
alcancen las condiciones de parada predeterminadas, 
proporcionando al final la mejor solución C* como 
salida. El pseudocódigo de la BVNS se describe en 
el Algoritmo 3. 


Después de la agitación se aplica la búsqueda lo¬ 
cal con la solución C como solución inicial. Esta 
búsqueda intenta eliminar cada una de las etique¬ 
tas de C', y luego agregar nuevas etiquetas siguien¬ 
do el método voraz MVCA, para reducir al míni¬ 
mo el número de componentes conexas, hasta que 


A. Ajuste de la amplitud de agitación 

El parámetro q m ax determina la amplitud que lle¬ 
ga a alcanzar el proceso de agitación en la VNS 
básica. El valor de este parámetro es relevante para 
modular la VNS básica y así conseguir un balance 
adecuado entre la potencia de intensificación y la ca- 















Algorithm 3: BVNS para el fcLSF 
Input: Un grafo no dirigido G = (V, A) 
con \V\ = n vértices y \A\ = m 
aristas que están etiquetadas con 
\L\ = í etiquetas y un entero po¬ 
sitivo k ; 

Output: Un bosque F; 

begin 

C* <— generar — al — azarQ ; 

repeat 

tomar q «- 1 y q m 

ax 

while q < q max do 

C <r- agitación(JV g (C*)); 

C" búsqueda-local(C"); 
if \C "| < \C*\ then 
moverse: C* «— C"; 

reinicializar la amplitud de agitación q <— 

i; 

actualizar la máxima amplitud de agita¬ 
ción qmax ^ \G |? 

else 

aumentar la amplitud de agitación: q <— 
9 + i; 

end 

end 

until condiciones de parada ; 

Tomar un bosque generador F del grafo 

G(C) = (V,A(C)) 

end 


pacidad de diversificación mi. La selección de un 
valor pequeño para q max produce una alta potencia 
intensificadora y una baja capacidad diversificadora, 
dando lugar a un algoritmo rápido, pero con una alta 
probabilidad de quedar atrapado en un mínimo local 
no óptimo. Por el contrario, un valor grande de q max 
disminuye la capacidad de intensificación y aumenta 
la capacidad de la diversificación, lo que da lugar a 
un algoritmo más lento, pero con mayor facilidad pa¬ 
ra escapar de los mínimos locales subóptimos. Con 
objeto de estudiar la relevancia del parámetro y op¬ 
tar por la mejor alternativa se establecieron tres ti¬ 
pos de estrategia. 

1. La primera estrategia consiste en elegir el paráme¬ 
tro qmax fijo para todas las ejecuciones del algoritmo 
( q m ax = a), independientemente de la instancia a la 
que se aplica. Por tanto, la potencia exploradora de 
la agitación es fija. 

2. La segunda estrategia consiste en elegir el 
parámetro q m ax dependiendo de la dificultad o mag¬ 
nitud de la instancia. En este caso optamos por to¬ 
mar la cota superior k del número de etiquetas a usar 
como una medida de la que depende directamente el 
tamaño del espacio de soluciones. Por tanto aplica¬ 
mos una fórmula del tipo q m ax = a ■ k. La potencia 
exploradora de la agitación se hace depender del ta¬ 
maño del espacio de soluciones. 


3. La tercera estrategia consiste en elegir el paráme¬ 
tro q m ax dependiendo del estado en que se encuen¬ 
tre la búsqueda. Como reflejo del estado en que se 
encuentra la búsqueda elegimos el tamaño de la so¬ 
lución actual. La fórmula para establecer el tamaño 
de la agitación es q max = a • |C|, siendo C la solu¬ 
ción actual. Por tanto, la potencia exploradora de la 
agitación se hace depender del estado en el que se 
encuentra el proceso de búsqueda medido en el valor 
de la función objetivo. 

Para evaluar el rendimiento de las estrategias des¬ 
critas para la BVNS se utilizó un subconjunto de 
las instancias propuestas en la literatura para este 
tipo de problemas (ver bd- Se trata de conjuntos 
de instancias generadas aleatoriamente a partir del 
número de nodos (n) y el número de etiquetas (£). 
Se obtuvieron conjuntos de instancias para diferentes 
combinaciones de los parámetros n, d y £. El número 
n de vértices varia de 100 a 500 y el número í de eti¬ 
quetas de 0, 25n a l,25ra. Para establecer el número 
de aristas se fijó la densidad del grafo en d = 0,5 pol¬ 
lo que el número máximo de aristas de cada instan¬ 
cia es m = 0,5n(n —1). El valor del parámetro k para 
adaptar estas instancias al fcLSFP se determinó ex- 
perimentalmente como \n¡ 2 J J, donde j es el valor 
más pequeño de tal manera que las instancias ge¬ 
neradas no formaran una única solución conexa una 
vez resuelto el MVCA [2¡. Para cada combinación 
de n, d, í y k se generaron al azar 10 instancias dis¬ 
tintas. La cota k para el número de etiquetas es un 
factor de importancia crucial para crear instancias 
significativas ya que, si se elige demasiado grande (en 
particular, mayor que el número de etiquetas necesa¬ 
rias para obtener un MLST), podría ser fácil para la 
heurística para detectar una sola componente conexa 
como solución óptima en un tiempo extremadamen¬ 
te corto. Por otro lado, si el valor de k es demasiado 
pequeño, el problema podría ser igualmente fácil, ya 
que el espacio de búsqueda podría ser recorrido muy 
rápidamente. 

Para analizar las etrategias diferentes de selección 
de la máxima amplitud de la agitación q m ax des¬ 
critas y optar por la que tenga mejor rendimiento 
se realizaron experimentos con 8 de estos conjuntos 
de instancias diferentes para cada una de las com¬ 
binaciones de los números de vértices n = 100 y 
n = 200 y números de etiquetas l = 0,25n, £ = 0,5n, 
£ = n, y £ = l,25n para un total de 80 casos. En to¬ 
dos estos casos se ejecutaron los algoritmos por un 
máximo de 60 segundos de tiempo CPU y el tiempo 
que aperece en las tablas corresponde al tiempo que 
emplearon en encontrar la mejor solución expresado 
en milisegundos. Las instancias están disponibles en 
http://www.sergioconsoli.com/MLSTP.htm. To¬ 
dos los cálculos se hicieron en un Mac OSX 10.9 con 
microprocesador a 3,0 GHz y con 8 GB de RAM. 

Se compararon 5 casos particulares de las tres 
estrategias anteriores eligiendo 5 valores para el 
parámetro a en cada una de ellas. Así para la 
estrategia 1 se consideraron los valores de a en 










{5,10,15, 20, 25}. Para la estrategia 2 se considera¬ 
ron los valores de a en {0,1,0,3, 0,5,0,7, 0,9}. Final¬ 
mente, para la estrategia 3 se consideraron los valo¬ 
res de a en {1/3,2/3,3/3,4/3,5/3}. Las 15 estrate¬ 
gias de actualización del parámetro q m ax correspon¬ 
dientes que dan lugar a 15 versiones de la BVNS se 
resumen en la Tabla [U Los resultados computacio- 
nales obtenidos con estas 15 versiones de la BVNS 
se muestran en las Tablas [Ti] a IIVI 


Num. 

Qmax 

Valores de a 

1 

= a 

5 

10 

15 

20 

25 

2 

= a ■ k 

0.1 

0.3 

0.5 

0.7 

0.9 

3 

= <*-\C\ 

1/3 

2/3 

3/3 

4/3 

5/3 


TABLA I: Estrategias de selección de la amplitud de 
agitación 

Los resultados computacionales obtenidos con ca¬ 
da una de las tres estrategia propuestas se pueden 
ver, respefctivamente, en la Tabla [Til en la Tabla. [TTT1 
y en Tabla ITVl 

Los resultados obtenidos con la estrategia 1 son 
bastante pobres aunque los valores del objetivo pa¬ 
ra los tamaños 20 y 25 son aceptables, aunque no 
el tiempo de ejecución. Los resultados Los resulta¬ 
dos computacionales obtenidos con la estrategia 2 
ya pasan a ser mejores pero se detecta una excesiva 
diversificación para valores grandes de a y los tiem¬ 
pos de ejecución son muy altos. Para la estrategia 
3 se observa un doble comportamiento. Cuando a 
es bajo, el algoritmo termina de manera prematura 
con soluciones de baja calidad, adoleciendo de capa¬ 
cidad de exploración. Cuando a es mayor que uno 
se observan mejores resultados porque la agitación 
se amplía. Sin embargo cuando a crece demasiado 
se produce un pico de diversificación dando lugar a 
malos resultados que se traducen en un tiempo ex¬ 
cesivo sin aportar soluciones de muy alta calidad. 
Como consecuencia de ello, tenemos que la tercera 
estrategia con una agitación de amplitud dependien¬ 
te del tamaño de la solución actual con una constante 
de proporcionalidad directa de 4/3, q m ax = 4|C|/3, 
produce los mejores resultados. Esta conclusión coin¬ 
cide con lo ya observado en ¡7¡ donde se muestra que 
el parámetro q m ax = (4|C|/3) da un buen balance 
entre exploración y explotación para el MLSTP. 

IV. Comparación entre metaheurísticas 

Para evaluar el rendimiento y eficiencia del al¬ 
goritmo propuesto en comparación con los descri¬ 
tos en la Sección m se usaron también las instan¬ 
cias propuestas en [2- Se realizaron experimentos 
con 20 conjuntos de instancias diferentes con un 
número de vértices n £ {100,200,300,400,500}, 
y número de etiquetas i £ {0, 25n, 0,5n, n, l,25n}. 
Para cada combinación de n y í, se dispone de 
10 instancias distintas, para un total de 200 ca¬ 
sos. Estas instancias también están disponibles en 


http://www.sergioconsoli.com/MLSTP.htm. 

Los resultados obtenidos se ofrecen en la Tabla [V] 
Las tres primeras columnas muestran los parámetros 
que caracterizan los diferentes conjuntos de instan¬ 
cias (n, ell y k ), mientras que las columnas res¬ 
tantes dan los resultados obtenidos con los algorit¬ 
mos descritos. Para cada conjunto de instancias, la 
calidad de la solución obtenida (Obj) se evalúa co¬ 
mo el número medio de componentes conexas entre 
las 10 instancias del problema. Un tiempo de CPU 
máximo permitido de 60 o 600 segundos fue elegido 
como condición de parada de todas los metaheurísti¬ 
cas, determinado experimentalmente con respecto a 
la dimensión de los conjuntos de datos considerados 
(es decir, 1 minuto para n < 200, y 10 minutos pa¬ 
ra n > 200). Para el algoritmo genético fijamos el 
número de generaciones para cada instancia de tal 
manera que los cómputos lleven aproximadamente 
al máximo tiempo de CPU permitido. Todas las me¬ 
taheurísticas se ejecutan hasta el tiempo máximo de 
CPU que se alcanza y, en cada caso, se registra la me¬ 
jor solución. Los tiempos computacionales de la Ta- 
bla[V]de las columnas etiquetadas con “Tiempo”son 
los tiempos medios en los que se obtienen las mejo¬ 
res soluciones. En el caso de que no se haya obtenido 
ninguna solución en 3 horas por el método exacto se 
denota por NF en la Tabla fVl El rendimiento de un 
algoritmo se puede considerar mejor que otro si ob¬ 
tiene menor valor objetivo promedio o con un valor 
objetivo promedio similar pero con un menor tiempo 
computacional. 

En la Tabla |V] se observa que la BVNS es capaz 
de obtener muy alto rendimiento para el problema 
kLSF. En todos los conjuntos de datos se obtiene 
muy rápidamente las soluciones con la mejor cali¬ 
dad. Con respecto al tiempo de ejecución, unas veces 
fue más lento que las otras heurísticas pero debido 
a que obtiene mucho mejores soluciones de calidad 
(por ejemplo, ver las instancias de [n = 200, í = 200], 
[n = 300, l = 375], [n = 400, í = 500], [n = 500, í = 
500]). Cuando la calidad media de la solución obte¬ 
nida con la BVNS fue comparable a la obtenida por 
los otros procedimientos, empleó un tiempo compu¬ 
tacional menor. En cuanto a las otras heurísticas, 
el algoritmo genético no tuvo muy buen rendimien¬ 
to tanto con respecto a la calidad de la solución y 
como con respecto al tiempo de cómputo, lo que con¬ 
firma, en general, un mal comportamiento de las me¬ 
taheurísticas poblacionales al abordar a esta clase de 
problemas (véase, por ejemplo 0, 0, i)- El Méto¬ 
do Piloto registró un mejor comportamiento que el 
algoritmo genético con respecto a la calidad media 
de las soluciones, pero los tiempos de ejecución eran 
aún muy grandes, lo que indica una dispersión exce¬ 
siva con respecto a la diversificación. GRASP parece 
ser una buena opción en cuanto al tiempo de ejecu¬ 
ción, pero la calidad media de las soluciones no fue 
tan buena como con la BVNS. Es interesante notar 
también que en todos los casos de problemas para 















Parámetros 

VNS-01-5 

VNS-01-10 

VNS-01-15 

VNS-01-20 

VNS-01-25 

n 

t 

k 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

100 

25 

3 

2.1 

20.5 

2.1 

1.2 

2.1 

0.7 

2.1 

0.9 

2.1 

1.9 

100 

50 

3 

5.2 

7.8 

5.2 

4.6 

5.2 

4.5 

5.2 

4.4 

5.2 

5.7 

100 

100 

4 

8.3 

49.8 

8.3 

24.6 

8.3 

164.7 

8.3 

19.9 

8.3 

85.8 

100 

125 

5 

4.1 

106.4 

4.1 

196.7 

4.1 

220.5 

4.1 

68.1 

4.1 

195.7 

200 

50 

3 

2.0 

128.6 

2.0 

37.9 

2.0 

109.6 

2.0 

98.9 

2.0 

45.7 

Eü 

100 

4 

2.2 

741.6 

2.2 

942.7 

2.2 

1085.4 

2.2 

797.5 

2.2 

933.2 

200 

200 

5 

7.0 

5379.6 

6.8 

4923.4 

6.8 

2072.8 

6.8 

4700.5 

6.8 

3580.9 

200 

250 

6 

4.4 

1218.3 

3.9 

9701.9 

4.2 

12522.3 

3.9 

10862.4 

4.0 

14854.4 


TABLA II: Resultados para la estrategia 1 


Parámetros 

VNS-02-1 

VNS-02-3 

VNS-02-5 

VNS-02-7 

VNS-02-9 

n 

£ 

k 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

100 

25 

3 

2.1 

1.5 

2.1 

1.4 

2.1 

1.0 

2.1 

1.1 

2.1 

1.6 

100 

50 

3 

5.2 

30.2 

5.2 

4.1 

5.2 

4.4 

5.2 

62.3 

5.2 

4.7 

100 

100 

4 

8.3 

68.9 

8.3 

125.9 

8.3 

44.1 

8.3 

73.8 

8.3 

83.2 

100 

125 

5 

4.1 

42.4 

4.1 

68.3 

4.1 

155.9 

4.1 

219.1 

4.1 

152.7 

200 

50 

3 

2.0 

118.9 

2.0 

69.7 

2.0 

38.2 

2.0 

127.9 

2.0 

17.4 

Ba 

100 

4 

2.2 

1943.9 

2.2 

1305.3 

2.2 

1407.6 

2.2 

1269.1 

2.2 

1068.9 

200 

200 

5 

6.8 

8298.7 

6.8 

3231.1 

6.8 

2156.3 

6.8 

3027.8 

6.8 

6570.6 

200 

250 

6 

4.1 

8638.9 

4.1 

11113.5 

4.1 

11183.2 

4.0 

17258.8 

4.0 

13368.5 


TABLA III: Resultados para la estrategia 2 


Parámetros 

VNS-03-1/3 

VNS-03-2/3 

VNS-03-3/3 

VNS-03-4/3 

VNS-03-5/3 

n 

e 

k 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

100 

25 

3 

7.5 

0.6 

7.5 

0.7 

2.1 

3.1 

2.1 

1.0 

2.3 

1.7 

100 

50 

3 

15.4 

0.6 

15.4 

0.8 

5.3 

2.6 

5.2 

4.2 

5.8 

4.5 

100 

100 

4 

29.3 

0.4 

00 

bo 

40.1 

8.7 

29.9 

8.3 

11.9 

8.3 

44.7 

100 

125 

5 

29.2 

0.7 

4.9 

572.4 

4.5 

23.9 

4.1 

36.1 

4.1 

38.8 

200 

50 

3 

9.0 

0.7 

9.0 

0.7 

2.1 

64.1 

2.0 

10.2 

2.3 

19.7 

m 

100 

4 

14.3 

0.2 

3.5 

362.4 

2.6 

658.2 

2.2 

731.1 

2.2 

2199.1 

200 

200 

5 

33.7 

0.5 

9.1 

238.1 

7.2 

673.9 

6.8 

1519.1 

6.8 

8636.5 

200 

250 

6 

36.2 

0.9 

6.3 

357.2 

4.6 

7979.3 

3.9 

2789.2 

4.1 

4600.1 


TABLA IV: Resultados para la estrategia 3 


los que el método exacto obtuvo la solución óptima, 
sólo la BVNS también produjo siempre la misma so¬ 
lución, y con un tiempo de cómputo muy corto. La 
BVNS proporciona el mejor compromiso entre las 
capacidades de intensificación y diversificación para 
el problema considerado. Resumiendo los resultados 
demuestran que la BVNS es una nretaheurística muy 
eficaz para el problema fcLSF, produciendo solucio¬ 
nes de alta calidad en tiempos computacionales de 
ejecución cortos. 

V. Conclusiones 

En este trabajo hemos considerado el problema 
del bosque generador /c-etiquetado (fcLSFP), un pro¬ 
blema NP-dwco que resulta de una extensión del 
problema del árbol generador de etiquetado míni¬ 
mo (MLSTP) al caso en el que se impone una cota 
superior k al número de etiquetas a utilizar; por lo 
que el propósito es obtener la solución con el mínimo 
número de componentes conexas que no violen esta 


restricción. Ya que el problema fcLSF es iV\P-duro, 
son interesantes las heurísticas y nretaheurística que 
permitan abordar con garantías de rendimiento. Se 
describe una BVNS para la que se ajusta el paráme¬ 
tro clave que fija la máxima amplitud que puede al¬ 
canzar el proceso de agitación. Se comprueba experi¬ 
mentalmente que la mejor opción es hacer depender 
dicho parámetro del estado del proceso de búsqueda, 
conclusión que podría extenderse a otros problemas. 
Al menos se presume que una fórmula de actuali¬ 
zación similar sería aplicable a los problemas en los 
que el tamaño de la solución refleja el estado del 
proceso de búsqueda. Con este ajuste del paráme¬ 
tro, la BVNS se comparó experimentalmente sobre 
una amplia gama de instancias del fcLSFP, con otras 
metaheurísticas propuestas en la literatura para el 
problema: el Método Piloto, un Algoritmo Genéti¬ 
co, y un procedimiento de búsqueda voraz aleatoria 
adaptativa (GRASP), junto con un enfoque exacto 
de resolución. Nuestros experimentos computaciona- 



















































































EXACTO 


PM 


GA 

GRASP 

BVNS 



k 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

Obj 

Tiempo 

100 

25 

3 

2,1 

59,4 

2,1 

6,8 

2,1 

47,2 

2,1 

7,7 

2,1 

1,1 

100 

50 

3 

5,2 

84,3 

5,2 

35,7 

5,2 

62,6 

5,2 

19,4 

5,2 

4,2 

100 

100 

4 

8,3 

1100 

8,3 

125,7 

8,3 

215,2 

8,3 

108,5 

8,3 

16,3 

100 

125 

5 

4,1 

44400 

4,3 

217,2 

4,1 

281,7 

4,3 

143,6 

4,1 

29,5 


50 

3 

2 

481,4 

2 

65,1 

2 

458,9 

2 

30,3 

2 

8,4 



4 

2,2 

2500 

2,4 

699,3 

2,7 

1100 

2,4 

319,1 

2,2 

243,6 

200 


5 

6,6 

978900 

6,9 

ywm 

7 

11300 

7 

775,2 

6,6 

1900 

200 

250 

6 

NF 

NF 

4,3 

MU 

4,4 

12700 

4,3 


3,9 

7500 

wm 

75 

3 

2 

793,8 

2 

295,5 

2 


2 

132,3 

2 

60,2 

m 

150 

4 

3,6 

16700 

3,8 

2600 

4 

IH 

3,9 

10100 

3,6 

4400 

300 

300 

5 

NF 

NF 

13,5 

12200 

13,6 

14100 

13,5 

13100 

12,9 

9700 

300 

375 

6 

NF 

NF 

10,3 

18900 

10,7 

16500 

10,1 

5900 

9,7 

41100 

400 

100 

3 

2,3 

1500 

2,3 

955,2 

2,4 

7200 

2,3 

365,1 

2,3 

236,9 

400 


3 

31,9 


31,9 

6500 

31,9 

10300 

31,9 

3600 

31,9 

503,8 

400 


5 

NF 

NF 

21,1 

38300 

21,6 

42900 

21,5 

7900 

20,9 

175000 

400 


6 

NF 

NF 

18,8 

58800 

18,7 

26900 

18,3 

33300 

17,7 

39500 


125 

3 

2,9 

21200 

2,9 

2200 

3,1 

22500 

2,9 

693,8 

2,9 

617,1 


250 

4 

5,8 

193900 

6,1 

19700 

6,3 

36800 

6,1 

9300 

5,8 


500 

500 

6 

NF 

NF 

8,8 


9,1 


8,8 

50300 

8,1 

166400 

500 

625 

7 

NF 

NF 

7,7 

155100 

8,9 


7,7 

140600 

7 

91700 

TOTAL: 

- 

- 

164,7 

426900 

168,1 

468800 

164,6 

287100 

159,2 

406600 


TABLA V: Resultados para |V| = 100,200,300,400,500 y \L\ = 0,25|V|, 0,5| V|, \V\, 1,25| V| 


les muestran que la BVNS debidamente ajustada su¬ 
pera a todas las otras metaheurísticas y es rápido, 
simple, y particularmente eficaz para el problema 
fcLSF. 

Agradecimientos 

Este trabajo está parcialmente financiado por el 
proyecto TIN2012-32608 del Ministerio de Economía 
y Competitividad. 

Referencias 

[1] T. Brüggemann, J. Monnot, G. J. Woeginger, Local 
search for the mínimum label spanning tree problem with 
bounded colour classes, Operations Research Letters 31 
(2003) 195-201. 

[2] Cerulli, R., A. Fink, M. Gentili and A. Raiconi, The k- 
labeled spanning forest problem, Procedía - Social and 
Behavioral Sciences 108 (2014), pp. 153-163. 

[3] Cerulli, R., A. Fink, M. Gentili and S. Vofí, Extensions 
of the mínimum labelling spanning tree problem , Jour¬ 
nal of Telecommunications and Information Technology 
4 (2006), pp. 39-45. 

[4] Cerulli, R., A. Fink, M. Gentili and S. Vofi, Metaheu- 
ristics comparison for the mínimum labelling spanning 
tree problem , in: B. L. Golden, S. Raghavan, E. A. Wa- 
sil (Eds.), The Next Wave on Computing, Optimization, 
and Decisión Technologies , Springer-Verlag, New York, 
2005, pp. 93-106. 

[5] Chang, R. S. and S. J. Leu, The mínimum labelling span- 
ninq trees, Information Processing Letters 63 (1997), 
pp. 277-282. 

[6] A. M. Chwatal, G. R. Raidl, Solving the mínimum label 
spanning tree problem by ant colony optimization , Pro- 
ceedings of the 7th International Conference on Genetic 
and Evolutionary Methods (GEM 2010) (2010) 91-97. 
Las Vegas, Nevada. 

[7] Consoli, S., K. Darby-Dowman, N. Mladenovic and J. A. 
Moreno-Pérez, Greedy randomized adaptive search and 
variable neighbourhood search for the mínimum labelling 


spanning tree problem , European Journal of Operational 
Research 196 (2009), pp. 440-449. 

[8] Consoli, S., K. Darby-Dowman, N. Mladenovic and J. A. 
Moreno-Pérez, Variable neighbourhood search for the mí¬ 
nimum labelling Steiner tree problem , Annals of Opera¬ 
tions Research 172 (2009), pp. 71-96. 

[9] S. Consoli, J. A. Moreno-Pérez, Solving the mínimum 
labelling spanning tree problem using hybrid local search, 
Electronic Notes in Discrete Mathematics 39 (2012) pp. 
75-82. 

[10] Goldberg, D. E., K. Deb and B. Korb, Don’t worry, be 
messy, in: Proceedings of the J^th International Confe¬ 
rence on Genetic Algorithms (1991), pp. 24-30. 

[11] Hansen, P. and N. Mladenovic, Variable neighbourhood 
search, Computers and Operations Research 24 (1997), 
pp. 1097-1100. 

[12] Hansen, P., N. Mladenovié, and J.A. Moreno-Pérez, Va¬ 
riable neighbourhood search: methods and applications, 
Annals of Operations Research 175 (2010), pp. 367-407. 

[13] Krumke, S. O. and H. C. Wirth, On the mínimum label 
spanning tree problem, Information Processing Letters 
66 (1998), pp. 81-85. 

[14] Miller, J. S., “Multimodal statewide transportation plan- 
ning: A survey of state practices,” Transportation Re¬ 
search Board, National Research Council, Virginia, US, 
2006. 

[15] Resende, M. G. C. and C. C. Ribeiro, Greedy randomi¬ 
zed adaptive search procedure, in: F. Glover and G. Ko- 
chenberger, editors, Handbook of metaheuristics, Kluwer 
Academic Publishers, Norwell, MA, 2003 pp. 219-249. 

[16] Xiong, Y., B. Golden and E. Wasil, A one-parameter ge¬ 
netic algorithm for the mínimum labelling spanning tree 
problem, IEEE Transactions on Evolutionary Compu- 
tation 9 (2005), pp. 55—60. 

[17] Xiong, Y., B. Golden and E. Wasil, Improved heuristics 
for the mínimum labelling spanning tree problem, IEEE 
Transactions on Evolutionary Computation 10 (2006), 
pp. 700-703. 

[18] Van-Nes, R., Design of multimodal transport networks: 
A hierarchical approach, Delft University Press, 2002 

















































